Magnetic properties of substitutional Mn in (110) GaAs surface and subsurface layers. 
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Motivated by recent STM experiments, we present a theoretical study of the electronic and mag- 
netic properties of the Mn-induced acceptor level obtained by substituting a single Ga atom in the 
(110) surface layer of GaAs or in one of the atoms layers below the surface. We employ a kinetic- 
exchange tight-binding model in which the relaxation of the (110) surface is taken into account. 
The acceptor wave function is strongly anisotropic in space and its detailed features depend on the 
depth of the sublayer in which the Mn atom is located. The local-density-of-states (LDOS) on the 
(110) surface associated with the acceptor level is more sensitive to the direction of the Mn magnetic 
moment when the Mn atom is located further below the surface. We show that the total magnetic 
anisotropy energy of the system is due almost entirely to the dependence of the acceptor level energy 
on Mn spin orientation, and that this quantity is strongly dependent on the depth of the Mn atom. 



I. INTRODUCTION 

Cross sectional scanning tunneling microscopy (STM) 
is a sophisticated and extremely valuable nanoscience ex- 
perimental tool. It enables, in particular, the manipula- 
tion and characterization of individual impurities in semi- 
conductors and metals with unprecedented spatial reso- 
lution and electronic sensitivity. In the last few years this 
technique has been used to study the electronic and mag- 
netic properties of individual Mn atoms and Mn atom 
pairs 2 ^ in GaAs. Interest in these studies stems in part 
from the notion that investigating Mn dopants in GaAs 
at the atomic scale could lead to a better understanding 
of and better control over magnetism in diluted magnetic 
semiconductors (DMSs). Recently STM methods were 
used to substitute individual Ga atoms with Mn atoms 
in the first layer of a GaAs(llO) surface. 3 4 STM allowed 
experimentalists to visualize the electronic properties of 
the acceptor wavefunction bound to the Mn position, and 
to probe the valence-band-hole mediated exchange inter- 
actions between two Mn atomsP Apart from its possible 
relevance in probing the basic physics of DMSs, this ex- 
periment is a remarkable example of how STM techniques 
can now be employed to engineer novel nanomagnets with 
specifically designed quantum mechanical properties. 

In Refs. |3|4j it was shown that both the Mn-induced 
(acceptor) in-gap state and the exchange interaction be- 
tween Mn pairs are strongly anisotropic with respect 
to crystallographic orientation. A theoretical analysis 4 
based on a tight-binding model of Mn atoms in bulk 
GaAs^l reproduces some of these features qualitatively. 
Nevertheless, it seems clear that the proximity of the 
Mn atoms to the surface must have a significant impact 
on their properties. In fact, a more recent study^, in 
which individual Mn atoms were carefully positioned on 
the GaAs(llO) surface and layer-by- layer on the first few 
layers below, showed that the reduced symmetry at the 
surface strongly modifies the wave function of the im- 
purity. The acceptor wavefunction properties depend on 



the precise substitutional depth at which the impurity is 
located. This experimental conclusion is supported by a 
recent theoretical studyP 

In this paper we present a theoretical study of in- 
dividual Mn dopants substituting for Ga atoms in a 
GaAs(llO) surface, or in one of atomic layers below the 
surface. Our aim is to provide a systematic analysis of 
how the electronic and magnetic properties of the ac- 
ceptor wavefunction are modified by the presence of the 
surface as the Mn impurity is inserted into successively 
deeper layers, and compare with the limiting case of a 
Mn in bulk GaAs. In contrast to previous studies^, we 
focus on the spin-orbit induced dependence of the accep- 
tor wave function on the direction of the Mn magnetic 
moment orientation. We compare our surface-influenced 
results with the results of Ref.0 in which the dependence 
of the acceptor wavefunction on the magnetic moment di- 
rection was studied for a Mn atom in bulk GaAs. 

Our analysis is based on a microscopic tight-binding 
model, which accounts for the crucial relaxation of the 
GaAs (110) surface layer and for spin-orbit interactions 
in the valence band which play an essentiaP role in 
(Ga,Mn)As magnetism. We do not account explicitly 
for the Mn (i-orbitals, but account for d — p hybridiza- 
tion instead by adding an effective exchange interaction 
between the Mn moment and valence band orbitals on 
nearest-neighbor As sites. We also include other inter- 
action terms to account for the Coulomb repulsion of 
electrons by the Mn ion. 

Our calculations show that the acceptor wavefunction 
is in general strongly anisotropic in space; the detailed 
anisotropy features depend very strongly on the sub- 
layer in which the impurity is located, in agreement with 
experiment 7 and with previous calculations. 8 For a rea- 
sonable choice of the parameters of our model, we find 
that for a Mn located in the topmost layer or in the first 
subsurface layer, the acceptor state has a large binding 
energy and a strongly localized wavefunction with a very 
weak dependence on the Mn magnetic moment direction. 
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As the impurity is inserted deeper beneath the surface, 
the acceptor wavefunction becomes progressively more 
delocalized and its dependence on the Mn moment ori- 
entation increases significantly. In particular, we find 
that the LDOS feature on the GaAs(llO) surface due to 
the acceptor level, the quantity which is probed most di- 
rectly by STM, is noticeably different between the cases 
of a magnetic moment pointing along the easy and hard 
magnetic directions when the impurity is located a few 
monolayers below the surface. This prediction could be 
tested in STM experiments in which the direction of the 
magnetic moment is manipulated with an external mag- 
netic field. We also show that total magnetic anisotropy 
is related in a simple way to the magnetic anisotropy 
of the acceptor state and that the magnetic anisotropy 
landscape depends in a non trivial way on the sublayer 
in which the impurity is located. 

Our paper is organized as follows. In Sec. [XT] we give 
an explicit description of the model we use, explaining 
its tight-binding Hamiltonian and elaborating on other 
details necessary to understand our findings. In Se c. [111 
we present the results of the model, starting in Sec. |III A 



with the case of a single Mn substituted for Ga in bulk 
GaAs. We then proceed in Sec. |IIIB] to discuss a substi- 
tutional Mn in a (110) surface layer, highlighting emerg- 



ing features and relating these to experiment. Sec. Ill C 



is devoted to the study of the transition from surface to 
bulk, in which the Mn is placed in successively deeper 
sublayers below the surface. Finally, we present our con- 
clusions and discuss the implications of our results in 
Sec. [IVJ The main text of this paper provides a detailed 
description of the Mn impurity at different depths below 
the (110) surface layer of the host semi-conductor. Some 
readers may wish to begin by reading the summary and 
conclusion section, which gives a brief description of our 
main results, before exploring the main body of the pa- 
per. 



II. THEORY 

In a (Ga,Mn)As III-V dilute magnetic semiconductor 
(DMS), the most energetically stable position for the Mn 
is a Ga atom site. This property can be understood 
as following from the atomic electronic structure of Mn, 
[Ar]3d 5 4s 2 , with the As 2 electrons allowing for the for- 
mation of crystal bonds similar to those between the 
Ga ([Ar]3d 10 4sV) and the As ([Ar]3d 10 4s 2 p 3 ) atoms 
of the host crystal. Because the Mn is missing the Ap 
valence electron of Ga, it acts as an acceptor. The Mn 
ion repels electrons and attracts a weakly bound hole, 
forming a neutral stateP^ Mn both generates the local 
magnetic moments (via the 3d 5 half-filled d-shell) and 
acts as a supplier of potentially itinerant holes th at can 
mediate their coupling. Zener's kinetic-exchangj- 9 * 11 * 12 * 
or indirect-exchange interaction applies to systems like 
(Ga,Mn)As in which local moments formed by the mag- 
netic impurities are coupled via itinerant s- or p-band 



carriers. 

Our study is based on tight-binding model with a 
Hamiltonian 



H — i^band + He 



Hso + i^coul • 



(1) 



that includes a kinetic exchange interaction between 
the local-moment and the band electrons. The kinetic- 
exchange model is appropriate^ when the Mn-d to As- 
p hopping amplitudes (see below) are smaller than the 
energetic separation between the <i-orbitals and the top 
of the valence band. The band term in Eq. [l ] is given 
in terms of the Slater-Koster parameter d 13 | 1 ^ ~for bulk 
GaAsP 



and 



EEE 

ij fi/j,' cr 



(2) 



Here, i and j are atomic indices, and \i and a are orbital 
and spin indices, respectively. The t 1 ^ , are the Slater- 
Koster parameters that do not depend on spin. The only 
non-zero parameters are the on-site energies (i = j, /i = 
//) and the nearest neighbor hopping matrix elements for 
the 5, p x , p y and p z orbitals. 

Our model contains s and p electrons only, the d- 
electrons of the Mn enter the model in the form of a 
Mn local moment with spin S = 5/2 spin which we 
treat classically in this paper. This local magnetic mo- 
ment is formed by the five 3<i-electrons of Mn that in the 
tetrahedral host results in bonding and antibonding sp-d 
states in the form 16 of a triplet of ^-symmetry (3d xy - 
, 3d zx - and 3^^-like) and an occupied doublet of e g - 
symmetry (3d x 2_ y 2- and 3d z 2-like). The doublet couples 
only weakly to the host, and is split from the triplet by 
the tetrahedral crystal field. The triplet hybridizes with 
the connecting sp-orbitals and the weakly bound hole oc- 
cupies one of the three sp-d antibonding states at the top 
of the valence band, predominantly of As 4p-character. 9 
The hybridization of the occupied Mn 3<i-electrons with 
the nearest neighbor As 4p-electrons, cause the p- states 
at the top of the valence band with spin parallel to the 
Mn spin to move up in energy relative those that are an- 
tiparallel, which hybridize with high-energy unoccupied 
d-orbitals. The direct exchange interaction between holes 
at the top of the valence band and the Mn ^-electrons is 
weak, such that p-d hybridization dominates which re- 
sults in an antiferromagnetic coupling ] 17 * 18 * This physics 
is captured in the H exc term of 0, which induces an 
exchange field on the nearest-neighbor As p-electrons, 



— Jpd ^2 n ' 



(3) 



m n[m] 



where J p d = 1.5 eV is the approximate value of the 
exchange coupling constant inferred from Refs. [T9] 
and [20], and 



y ] a n7T(j T o-a f 



(4) 
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where f is the Pauli matrix vector. The first sum over 
m in ([3| runs over all Mn, and the second sum runs over 
all As that are nearest neighbors to Mn atom m, de- 
noted by n[m\. Because the exchange field in this model 
influences only p-electrons, the sum ir in Q only runs 
over the three p-orbitals. The valence band electronic 
structure depends on the classical Mn spin orientation 
Cl m (which we parametrized by polar and azimuthal <\> 
angles) through the scalar product S n • Cl m . 

We approximate the spin-orbit coupling Hamiltonian 
-^so by a local atomic one-body operator, in which the 
spin quantization axis is defined by (l m : 



(5) 
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where i is an atomic index and A^ denotes the renor- 
malized spin-orbit splitting 15 for which we use the values 



A,- 



e{Ga} 



58 meV, A, 



ie{As} 



140 meV, and A 7 - 



e{Mn} 



Ke{Ga}/2- The spin-orbit term causes the total energy 
of the system, as obtained by summing the energies of 
all occupied eigenstates, to depend on the magnetization 
direction parametrized by Cl m . Our procedure allows us 
to calculate the magnetic anisotropy landscape on the 
unit sphere and extract the magnetic anisotropy energy 
as £?anis = max — min E (0, (/>). The polar coor- 

dinate system is defined such that = corresponds to 
the [001] direction, and with = 7r/2, </> = and <\> = ir/2 
correspond to the [100] and the [010] directions respec- 
tively. 

The presence of a negatively charged Mn ion at- 
tracts holes and repels electrons. We represent the spin- 
independent part of the effective potential due to Mn 
substitution by a long-range repulsive Coulomb part and 
a Mn central cell correction term, 



2 t 



(6) 



The first term in (|6| represents the long-range part, 
which is reduced by the host material dielectric constant 
e r = 12. To account crudely for weaker dielectric screen- 
ing at the surface, the dielectric constant for a Mn on 
the surface is reduced to e r = 6 for the affected surface 
atoms. The correction term consists of on- and off-site 
parts, V corr = V on + V f[ which influence the Mn ion and 
its nearest neighbors respectively. The on-site Coulomb 
correction is estimated to 1.0 eV from the ionization en- 
ergy of Mn. The off-site Coulomb correction affects all 
the nearest-neighbor As surrounding the Mn ion and to- 
gether with i^exc (|3| reflects primarily p-d hybridization 
physics. It is one of the most important parameters of 
the model and its value is set by tuning the position of 
the Mn-induced acceptor l evel in the bulk to the exper- 
imentally observed position ^ 1 * 22 * 23 * 24 ! at 113 meV above 
the first valence band level. The value thus obtained is 
V fi = 2.4 eV. In this picture, long-range Coulomb, ex- 
change and correction interactions all play an important 
role in determining the character of Mn acceptor levels. 



FIG. 1: The relaxed (110) surface. Gray filled circles signify 
relaxed positions and the black filled circles the unrelaxed 
positions. This illustration shows a side view of the (110) 
surface with distances in units of the GaAs lattice constant 
aiatt — 0.565 nm. 
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FIG. 2: (Color online) The (110) surface. The green spheres 
represent Ga, blue spheres the As and red sphere indicates a 
Mn that has replaced a Ga. 



We model the electronic structure of GaAs with a sin- 
gle substitutional Mn by performing a super-cell type cal- 
culation with a cubic cluster of 3200 atoms and periodic 
boundary conditions in either 2 or 3 dimensions, depend- 
ing on whether we are studying the (110) surface or a 
bulk-like system. The (110) surface of GaAs is simplified 
from both theoretical and experimental points of view, by 
the absence of large surface reconstruction. Relaxation 
of surface layer positions must nevertheless be included 
since it removes dangling-bond states that would oth- 
erwise obscure the band-gap. We follow the procedure 
outlined in Refs. [25 26J, in which atomic shifts as deep 
as the 2nd sublayer are taken into account. The resulting 
(110) surface is summarized in Fig. T] and the top layer 
and first subsurface layer are depicted in Fig. [2] The most 
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salient feature is an up-shift of the surface As, accompa- 
nied by a down-shift of the surface Ga. Both species are 
shifted along the surface, such that the effective bond- 
angles change, which affects the direction cosines entering 
the Slater-Koster parametrization. 13 The tight-binding 
parameters scale like d^/d 2 , where do is the original dis- 
tance and d is the distance after the relaxation. The 
rescaling parameter is small and is at most d^/d 2 ~ 1.06. 

By projecting the eigenvector of the acceptor level onto 
surface sites obtained after diagonalizing the Hamilto- 
nian 0, we can study the surface LDOS, which is inti- 
mately related to the topographic STM images:^ In our 
cluster approach we sum the absolute square of the co- 
efficients in the eigenvector corresponding to the LDOS 
on a given atom j, 

(V-nl^hM = E ^^{i'fi'a'lPjlifia) (7) 

= EI c "m.| 2 > (8) 
no- 
where n is an eigenvalue index, i and j are atomic indices, 
and \± and a denote orbital and spin respectively. In ^ 
Pj = X^cr b/^XjA^ 7 ! projects out the LDOS of atom j. 
Similarly, we can define operators that project out the 
orbital and spin character for a given eigenlevel n: 

p » = K/^X^H > ( 9 ) 

The procedure we follow in generating LDOS plots is sim- 
ilar to the approach used by Tang et al. in Ref. |5] We 
place Gaussians on the atomic positions with a magni- 
tude equal to the LDOS of that atom, and a full-width 
at half-maximum equal to half the nearest neighbor spac- 
ing. (This procedure mimics the finite spatial resolution 
of an STM tip). In all the LDOS plots the normaliza- 
tion is such that the sum of LDOS over all atoms in the 
cluster is unity for a single eigenlevel. When the STM 
is operating in constant current mode (see for example 
Ref. |27j) , the tunneling current is maintained at a fixed 
value by varying the tip-surface distance. The exponen- 
tial decay of the surface wavefunction causes the distance 
recorded to depend approximately logarithmically on the 
LDOS at the surface. We therefore employ a logarithmic 
color-scale in our images of the LDOS. 

When imaging states in a semi-conductor band gap 
at a low bias, it is necessary to move the tip very close 
to the substrate surface. This means that interactions 
between the tip and the sample can cause a change in 
the surface wavefunction. In addition, many-body ef- 
fects beyond those captured by the mean-field descrip- 
tion of electronic states outlined above can be important 
in some cases. In particular, as detailed below, we find 
that acceptor levels for Mn very close to the surface lie 
deep in the gap. STM experiments provide a partial pro- 




FIG. 3: (Color online) In (a) the magnetic anisotropy energy 
of a single Mn in a 3200 atom cluster with periodic boundary 
conditions. Bistable easy directions parallel to the [001] axis 
are separated by a single barrier of magnitude 4.35 meV in the 
(001) plane, (b) shows the ground-state energy level spectrum 
in the easy direction, where the off-site Coulomb correction 
term been used to tune the acceptor level at 113 meV above 
the first valence band level. 

file of the spatial distribution of added or removed elec- 
trons. When the acceptor is deep in the gap, addition 
or removal will cause a big change in the potential seen 
by other electrons. These various many-body effects are 
partly captured by our phenomenological model, but we 
cannot expect to find exact correspondence with the ex- 
perimental images. 



III. RESULTS 
A. Single Mn in bulk GaAs 

Our starting point is a single Mn in bulk, as repre- 
sented in our model by placing it at the center of a cu- 
bic GaAs cluster of 3200 atoms and enforcing periodic 
boundary conditions. The dimensions of the cluster (in 
the crystalline directions [110] x [110] x [001]) used in all 
calculations are 38.0 x 38.0 x 42. 4A 3 , or in terms of atomic 
layers 20 x 20 x 32. Creating the supercell in the form 
of a cubic cluster defined by the limiting planes (110), 
(110) and (001), enables us to apply periodic boundary 
conditions in two directions and study the (110) surface. 
The generated cluster can also be given periodic bound- 
ary conditions in all three directions to create a 'bulk' 
system, which should simply be regarded as the fully pe- 
riodic counterpart of the surface system. 

The magnetic anisotropy energy for the fully peri- 
odic system as a function of the magnetization direction 
E(Q, </>), is shown in Fig. [3] We find bistable minima and 
an easy axis parallel to the [001] direction, separated by a 
single barrier equal to £^ a nis = 4.35 meV. The anisotropy 
is very sensitive to the cubic symmetry and the particular 
easy direction [001] can be seen as a consequence of the 
supercell symmetry. The closest distance between Mn in 
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adjacent supercells is 38 A along the [110] and [110] direc- 
tions, along which the most prominent hopping path oc- 
curs via closely spaced As and Ga. The distance along the 
4 other equivalent symmetry directions is much longer, 
resulting in the observed anisotropy. Scaling down the 
size of the supercell (and thereby increasing the effective 
Mn doping) to 1200 atoms yields an anisotropy of 8.8 
meV, and the value increases further as the supercell size 
is decreased. We conclude from this calculation that the 
magnetic anisotropy of a single Mn in the bulk of an in- 
finite crystal is much smaller than ~ 4 meV, and that, 
as far as magnetic anisotropy is concerned, the dilute 
isolated impurity limit is achieved only at Mn atom frac- 
tions x much smaller than 10 -3 , and much smaller than 
what can be represented in this or any other supercell 
calculation. 

Fig. [4] shows plots of the LDOS of the acceptor level in 
the easy and hard directions. Sim ilar res ults have been 
obtained previously by Tang et al. | 4 | 5 | 6 | T | using a similar 
tight-binding model. When the Mn spin is pointing in the 
easy direction [001], the wavefunction spreads out sym- 
metrically along [110] and [110]. The spread along [110] 
is shown is shown in Fig. El (a)-(e) in the left column, 
in which cuts of successive (110) planes up to 4 layers 
away from the Mn plane are displayed. These images 
agree qualitatively with the results presented in Ref. [TJ 
in which statistical methods and comparison with theory 
enabled the identification of the location of a particu- 
lar single Mn down to the fourth sublayer. The center 
and the right column in Fig. [4] show the LDOS in suc- 
cessive (110) planes when the Mn spin is oriented in the 
[110] and [110] hard directions respectively. The images 
show that the acceptor wavefunction exhibits a definite 
preference to spread out along the symmetry direction 
perpendicular to the Mn spin. In the center column of 
Fig. [4] it can be seen that the spread along the Mn spin 
direction [110] is very weak. By contrast, when the Mn 
spin is pointing in the [110] direction, perpendicular to 
the (110) plane normal, the rightmost column of Fig. [4] 
reveals much higher values of the LDOS. 

The cut of the Mn plane in the easy direction (Fig. [4] 
(a)) has 6% of the spectral weight of the acceptor wave- 
function on the Mn atom, a total of 15.4% on the 4 
nearest neighbors As and the rest is spread out in the 
lattice. These values correspond to a more spread out 
LDOS than the one obtained by Tang et al. in Ref. 
who find 10% on the Mn and 20% on the 4 surround- 
ing As. In both cases model parameters were adjusted 
to give the correct energetic position for the bulk Mn 
acceptor level. The difference in wavefunctions demon- 
strates that satisfying this criterion does not guaran- 
tee that the character of the acceptor level is correctly 
captured. In particular, accounting for the contribu- 
tion to binding from longer ranged Coulomb interac- 
tions lead to more extended wavefunctions at a given 
acceptor energy. The observed correlation between spa- 
tial anisotropy and Mn spin direction agrees qualitatively 
with previous calculations^, in which a decrease of 90% 
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FIG. 4: (Color online) The bulk acceptor level LDOS of a 
single Mn in a cluster with periodic boundary conditions in 
all directions. The figures show (110) plane cuts of the LDOS 
containing the Mn in (a), and moving from one to four lay- 
ers away from it in (b)-(e). The direction of the Mn spin 
is indicated by white arrows at the top of the figure. The 
left column shows the LDOS when the Mn spin is in the easy 
[001] direction, and the center and right column when the Mn 
spin is in two hard directions, [110] and [110]. The successive 
cuts reveal a decreasing magnitude of the LDOS with dis- 
tance from the impurity and a strong dependence on the Mn 
spin orientation. Note the alternating behavior on odd and 
even layers away from the Mn plane. When the Mn spin is 
in the easy direction, the acceptor wavefunction spreads sym- 
metrically in the [110] and [110] directions. Comparing the 
case when the spin is in a direction parallel to the (110) plane 
normal (center column), and when the spin is perpendicular 
to the (110) plane normal (right column), it can be seen that 
the wavefunction is highly anisotropic and extends in along a 
symmetry direction perpendicular to the spin direction. 



in maximum spectral weight at the center of the (110) 
plane images is seen four layers away from the Mn (cor- 
responding to Fig.[4](e)), when the Mn spin changes from 
[110] to [110], and a decrease of 15% when the Mn spin 
changes from [110] to [001]. Comparing this with our re- 
sults we find similar anisotropics. To begin with consider 
the anisotropy of the wavefunction at three layers away 
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FIG. 5: The total magnetic anisotropy energy obtained by 
summing the dependence of energy on magnetization orien- 
tation for all occupied orbitals as a function of the index of 
the Fermi level orbital. This figure shows that all orbitals in 
the valence band are dependent on magnetization orientation, 
causing the anisotropy to fluctuate substantially when one oc- 
cupation number is changed. The systems total anisotropy 
energy is indicated by a filled circle. We also see that the 
anisotropy is nearly zero when all valence band orbitals are 
occupied (indicated by empty circle). The anisotropy energy 
of a Mn acceptor is therefore just the anisotropy of the top 
level in the valence band, with a change of sign. 




FIG. 6: (Color online) The acceptor level and p-character 
variation on the unit sphere. The acceptor level anisotropy 
(a) cancels the anisotropy built up by all occupied levels. This 
means that the variation of the acceptor level as a function 
of the magnetization direction is the negative of the total 
anisotropy landscape. Variations in the orbital character (c-d) 
of the acceptor correlate with an acceptor wavefunction that 
spreads mainly in direction perpendicular to the Mn spin. 



from the Mn (Fig. [4] (d)). Here the maximum LDOS in 
the (110) plane drops by 86% as the spin changes from 
[110] to [110], and by 21% when the spin changes from 
[110] to [001]. Looking at the next layer (Fig. [5] (e)), 
we note that this anisotropic effect decreases slightly, 
and we find that the maximum LDOS decreases by 74% 
when the spin direction is changed from [110] to [110], 
and by 25% when it is changed from [110] to [001]. The 
source of this difference in the behavior of the LDOS be- 
tween odd and even layers away from the Mn, is that 
the odd layer maximum is on the As and the even layer 
maximum is on the Ga in the nearest neighbor hoping 
path along the [110]. This effect becomes pronounced 
three layers away from the Mn and higher. Although 
the actual percentages differ somewhat, the same strong 
anisotropic behavior of the acceptor wavefunction as in 
Ref. [6] is observed. The strong dependence of the LDOS 
on Mn spin-orientation, combined with the weak mag- 
netic anisotropy, implies that the observed LDOS should 
be strongly sensitive to external magnetic fields. 

The anisotropy energy is formed by summing up the 
energies of all occupied levels, counting 4 electrons per Ga 
and As and 3 electrons per Mn. The spin-orbit induced 
level shifts vary with magnetization direction on the unit 
sphere and give rise to an anisotropic dependence of the 
total energy on the magnetization direction. Fig. [5] shows 
the anisotropy energy as obtained by successively sum- 
ming eigenlevel anisotropy landscapes on the unit sphere, 
starting from the lowest level. It reveals fluctuations that 
persist far into the valence band. Fig. [5] illustrates the 
advantage of using a hole rather an electron picture in 
analyzing the anisotropy; the anisotropy built up by the 



shifts of many occupied levels, is retrieved by the single 
unoccupied acceptor level. This picture remains valid as 
long as the coupling to the conduction band is not sen- 
sitive to magnetization orientation. Quite generally we 
find that the anisotropy of the acceptor level e acc (0, (j)) 
and the systems total energy E(0, (j)) are accurately re- 
lated by 

e aC c (0, (j>) - Min[e acc (0, <j>)] 
= -(E(9^)-Min[E(9^)}). (11) 

This relation holds for several Mn in the system, subject 
to the same conditions. The single level e acc (0, (j)) should 
then be replaced with a sum over all acceptor levels. Be- 
cause the number of occupied levels that contribute to 
the total anisotropy varies from case to case, the hole 
picture is invariably more useful in trying to understand 
trends. 

Fig. [5] shows that the anisotropy energy built up by 
all occupied levels (filled circle), is canceled by adding 
the anisotropy of the single acceptor level (empty circle) . 
In Fig. [6] (a) the acceptor level variation as a function of 
magnetization direction is shown, and it can be seen that 
taking the negative of this reproduces the total energy 
variation shown in Figj3j The property observed in the 
LDOS images (see Fig. Hi) that the acceptor wavefunction 
tends to spread out in a plane perpendicular to the Mn 
spin direction is reflected in the orbital character varia- 
tion of the acceptor level. Fig.[6](b)-(c) shows the orbital 
p-character of the acceptor level, as obtained by using the 
projector To illustrate this point, consider the Mn 
spin in the direction [001], defined by = 0. Here, the 
p^-character of the acceptor wavefunction is lowest and 
the p x and p y characters are high, consistent with an ac- 
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ceptor wavefunction that is spreading out in the (001) 
plane. In the [100] direction, defined by 6 = tt/2 and 
= 0, the p x character is low and the p z and p y char- 
acters are high, as the wavefunction is spreading mainly 
in the (100) plane. A similar dip in the p^-character is 
found in the [010] direction. The sum of the p x and p y 
characters vary in an opposite manner to p zi such that 
Px +Py — 0.9 — Pz-> which means that p x +p y is constant 
along the line — tt/2. In the [110] hard direction where 
6 = 7r/2 and ip = 7r/4, the wavefunction should extend 
mainly in the (110) plane, which can then be seen by 
approximately equal characters of p x , p y and p z . 

Because of the spin-orbit interaction, the levels do not 
have definite spin. The highest occupied level acquires a 
minority-spin component that varies between 6-8% while 
the acceptor has a much smaller minority spin compo- 
nent between 0.2-0.3%. The shallower the impurity level, 
the more minority-spin character is acquired; the second 
highest occupied level has larger minority spin character 
in the range 12-16%. 



B. Single Mn in the (110) GaAs surface layer 

We now turn to the case of a single Mn in the (110) 
surface layer which has three As nearest neighbors, two 
located on the surface - see Fig. [2] In the calcula- 
tions we use the parameters in the Hamiltonian that 
were obtained by fixing the correct acceptor binding en- 
ergy in bulk. The magnetic anisotropy landscape on the 
unit sphere is shown in Fig. [7] (a). At the surface, the 
anisotropy energy range is E an i s = 3.5 meV and the land- 
scape has bistable minima and an easy axis at an ap- 
proximate 45° angle to the (110) surface, corresponding 
to the [111] direction. Fig. [7] (b) shows a portion of the 
eigenlevel spectrum. We find that the acceptor level of a 
surface Mn is very deep in the gap. The surface Mn ac- 
ceptor level is not at all similar to its bulk counterpart, in 
sharp contrast with what has previously been assumed. 
The highest occupied level is also deep in the gap 12-19 
meV below the acceptor level. The loss of coordination at 
the surface is primarily responsible for the much deeper 
state. A deep acceptor is also observed in experiment 4 , 
where the dl/dV curve reveals a broad resonance at 850 
meV above the valence band edge. Our results imply 
that the acceptor at the surface can be categorized as an 
intrinsically deep state, although the exact position de- 
pends on the degree of p-d hybridization at the surface. 

On the basis of generic considerations which recog- 
nize the reduced symmetry at the surface, one might 
have expected that the anisotropy would be highest at 
the surface. Due to the nature of the Mn interaction 
with neighboring As, we find that this is not the case. 
A high anisotropy requires a hole wavefunction that is 
more spread out in the lattice causing high variations in 
the orbital and spin character of the hole state. Fig. [7] 
(c)-(d) reveals variations in p x and p y orbital character 
of 1% and just a fraction of a percent in p z character. 




FIG. 7: (Color online) The magnetic anisotropy energy of a 
single Mn on the (110) surface layer. Panel (a) shows the 
magnetic anisotropy energy with bistable easy directions at 
an approximate 45 degree angle to the (110) surface in the 
[111] direction, separated by a barrier of 3.5 meV. In (b), 
on the left, the eigenlevel spectrum in the easy direction is 
shown, with the highest occupied level indicated by a red 
line. The surface acceptor level is now deep in the gap at 
1.27 eV. Reducing the off-site Coulomb correction to 1.57 eV, 
causes the acceptor level to lie less deep at the experimental 
position of the resonance at 0.85 eV (right hand spectrum 
in (b)). The anisotropy correlates with the variation in the 
acceptor orbital and spin character, shown in panels (c)-(f). 



This is in turn connected with the LDOS of the hole that 
is more localized than bulk. Fig. [8] (a) shows the ac- 
ceptor LDOS of the surface layer. The maximum spec- 
tral weight of 23.0% is not located on the Mn (which 
has 12.0% spectral weight), but rather on its 2 surface 
As nearest neighbors. The third As nearest neighbor in 
the sublayer has a much smaller spectral weight of 5.2%, 
which means that the core region of the Mn and its 3 
neighbors contains a total of 63% spectral weight. The 
Ga atoms have a smaller maximum spectral weight of 
6.1% for the top layer and 2.0% in the sublayer (see Fig. [8] 
(b)). Although the LDOS is highly localized, a similar- 
ity with the star-shaped symmetric images observed in 
the STM topographs 4 7 can be seen. The presence of 
Zn dopants in experimental samples might increase the 
coupling to conduction states, causing a more extended 
surface wavefunction. The surface and subsurface LDOS 
is similar in easy and hard magnetization directions; the 
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FIG. 8: (Color online) The LDOS for a Mn at the (110) sur- 
face, (a) shows the LDOS at the surface and (b) in the sub- 
surface when the Mn spin is pointing in the easy direction. 
The acceptor LDOS is much more localized than the bulk 
acceptor LDOS, with a total of 63% spectral weight on the 
Mn and its 3 neighboring As. (c) shows the effect of reduc- 
ing the off-site Coulomb correction in order to reproduce the 
surface acceptor level at the experimentally observed position 
of 850 meV. Although this change shifts energies appreciably, 
it produces only very small changes in the spatial pattern of 
the LDOS. The star-like shape found in experiment!^ can be 
distinguished in (a) and (c). (d) shows a maximal intensity 
projection (see main text for definition) of the whole cluster 
along the line of sight (the [110] direction) and reveals the 
symmetry responsible for the minima along [111]. 



changes in energy with magnetization direction are due 
to the changes in p x: p y and spin character on the same 
atoms. 

The position of the acceptor state is sensitive to the off- 
site Coulomb correction. Setting this parameter to repro- 
duce the experimental position gives an offsite Coulomb 
correction V q r = 1.57 eV. The qualitative properties 
of the state does not change but the magnitude of the 
anisotropy drops to 0.87 meV as we move away from the 
conduction band and the gap between the highest occu- 
pied and the acceptor level increases. Fig. 8 (c) reveals 
that this state has a very similar acceptor surface LDOS. 
Our calculations indicate that the acceptor state for a 
single Mn on the (110) surface is a deep, highly localized 
state with a relatively low anisotropy. Precisely how deep 
and how localized this acceptor level is, is dependent on 
model parameters that we must choose phenomenolog- 
ically. It is possible that the experimental surface Mn 
acceptor level binding energy quoted above could be in- 
accurate because of band-bending effects which could be 
present when performing STM on a semiconducting sur- 
face. Experimental studies of surface Mn in more heav- 
ily doped (Ga,Mn)As samples, in which band-bending 
effects are weaker could help settle this question. 

We can understand the occurrence of the bistable min- 
ima approximately along [111] from Fig. [8] (d), which 
shows a maximal intensity projection from a side view 
of the Mn, along the [110] direction. The presence of 



the surface, causes the wavefunction to spread out in the 
(111) plane down into the lattice. The star-like protru- 
sions across rows observed in experiment are weak sur- 
face echoes of this spread. The maximal intensity pro- 
jections are generated as follows. Instead of placing two 
dimensional Gaussians at the atomic sites (like in Fig. [8] 
(a)-(c)), we place a Gaussian sphere at each atomic site 
xf ^e^ i_r ) / r . Here, the magnitude of the Gaussian 
li is equal to the LDOS value of atom i (as obtained 
from Eq. [7]), and we choose the Gaussian smearing T 
such that the full width at half maximum is equal to half 
the nearest neighbor distance. This generates a three 
dimensional LDOS view of the entire supercell cluster, 
with one color value corresponding to the LDOS for each 
point in space. The maximal intensity projection takes 
the maximum LDOS value along the viewers line of sight 
(in this case the [110]) and projects it onto the viewing 
plane, i.e. the plane which is perpendicular to the view- 
ers line of sight. In this way, we obtain an image which 
is not directly related to the STM images, but provides 
useful information on the LDOS below the surface. 



C. Single Mn in subsurfaces layers of GaAs(llO) 

This subsection is devoted to the study of what hap- 
pens when the Mn is successively moved down from the 
surface layer, toward the cluster center layer which best 
approximates the bulk. Fig. [9] shows how the anisotropy 
energy varies over the unit sphere with increasing depth, 
with panels (a)-(i) corresponding to sublayers 1-9. The 
subsurface Mn landscape in Fig. [9] (a) reveals an ex- 
tremely low anisotropy energy of 30/ieV. As in the sur- 
face layer case, the small anisotropy can be traced to a 
very low variation in orbital character associated with a 
highly localized hole wavefunction. In the second sub- 
layer (Fig. [9] (b)), the Mn produces an anisotropy land- 
scape similar to that of the surface, but with high and 
low blocking barriers of 2.2-2.9 meV. In the successive 
layers, the high and the low barriers both increase and 
become larger than bulk, reaching a maximum of 14.5- 
10.9 meV in sublayer 5 (Fig. [9] (e)). In sublayer 6 (Fig. [9] 
(f)) the high and the low barrier have interchanged po- 
sitions. The easy axis remains approximately along the 
[111] all the way down to sublayer 6 where it begins to 
shift towards the surface normal. The low barrier de- 
creases towards the center and at the deepest layer of the 
slab the low barrier becomes so low that an approximate 
easy (110) plane forms (see Fig. [9] (i) at <j) = 7r/4,57r/4), 
opening up reversal paths connecting the two bistable 
minima. In the deepest sublayers the Mn spin can ex- 
plore the plane with <j> = 7r/4,57r/4 and any 6 with a 
very low energy cost. This plane corresponds to (110), 
which can then be classified as a quasi-easy plane. At the 
deepest levels, we still see traces of the now very shallow 
surface minima. 

(a) shows the evolution of the three highest 
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occupied levels and the acceptor level as a function of 
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FIG. 9: (Color online) Magnetic anisotropy over the unit 
sphere as a function of depth, (a)-(i) correspond to the 
anisotropy landscapes of sublayers 1-9. Sublayer 1 has ex- 
tremely low anisotropy, with easy directions along the [110]. 
In sublayer 2 (panel (b)) a surface-like landscape reappears, 
but with a high and a low blocking barrier. The low barri- 
ers are marked with an open circle, the high barrier (hard) 
directions with a filled circle and the minimum energy easy 
directions with a square. The barriers both grow with depth 
until they exchange positions in sublayer 6. The anisotropy 
energy reaches a maximum of 15 meV in sublayer 5 in (e). The 
difference between the high and the low barrier then increases 
with depth until a quasi-easy (110) plane forms (marked with 
dashed line), (j) tracks the high and the low barrier as a func- 
tion of the Mn depth, where the horizontal line indicates the 
bulk blocking barrier. 
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FIG. 10: Energies of the three highest occupied levels (filled 
circles) and the acceptor level (empty circles) as a function 
of Mn depth, starting from a Mn on the surface (layer 0), 
when the Mn spin is pointing in the easy direction (marked 
by squares in Fig. [9]). The acceptor abruptly drops down 
towards the valence band in sublayer 2 and then flattens out, 
converging at 125 meV above the first valence band level, 
corresponding to a slightly deeper acceptor than in bulk, (b- 
c) Energies of these four levels as a function of sublayer index, 
starting from sublayer 2; (b) is for the easy direction and (c) 
for the hard direction, respectively. The hard direction is the 
magnetization direction corresponding to the high barriers 
(marked by filled circles in Fig. |9j). 



Mn depth, when the Mn spin is pointing in the mag- 
netic easy direction (marked by squares in Fig. [9}. The 
acceptor abruptly drops down towards the valence band 
in sublayer 2 and then flattens out, converging at 125 
meV above the first valence band level, corresponding to 



a slightly deeper acceptor than in bulk. In Fig. [TO] close- 
ups of the same energy levels is shown in the easy (b) 
and hard (c) direction, respectively. The hard direction 
is given by the magnetization direction corresponding to 
the high barriers marked by a filled circle in Fig. [9j 

The difference between the acceptor levels in the 
hard and easy direction is the spin-orbit induced total 
anisotropy energy. The maximal anisotropy in sublayer 
5 corresponds to a very large variation in the gap be- 
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tween the highest occupied level and the acceptor of 5-50 
meV between the hard and the easy direction. The two 
quasi-degenerate states are split by the spin-orbit shift. 
The easy direction is the direction in the which the gap 
caused by this shift, is maximized, leading to a decrease 
of the total energy of the system. A plot of the high 
and the low barrier as a function of sublayer depth is 
shown in Fig.[9](j). From this figure we see the magnetic 
anisotropy energy is maximally enhanced in sublayer 5. 
This is due to the presence of the surface and has to do 
with the way the acceptor wavefunction extends around 
the Mn, as we move down through the layers. 
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Fig. 11 shows how the LDOS in the (110) surface plane 



evolves as a function of the sublayer depth in the hard 
and easy direction respectively. This figure therefore re- 
lates directly to the window on acceptor level proper- 
ties opened by STM experiments. The left column of 
Fig. 11 shows maximal intensity projections (see defi- 
nition above) for the easy direction solution. This se- 
ries of images gives a qualitative idea of where the main 
spectral weight of the acceptor wavefunction is located 
in relation to the surface (indicated by dashed line). We 
see that when the Mn impurity is close to the surface, 
the buckling has a large effect on its spatial extension. 
As the acceptor wavefunction detaches from the surface 
with increasing depth, the extension towards the surface 
becomes reduced and it begins to extend along [110] (not 
shown in this sequence of images as the extension along 
[110] is parallel to the line of sight). 

The center and the right column of Fig. [TT] shows the 
(110) surface LDOS when the Mn spin is pointing in the 
easy and hard direction, as depicted in Fig. [9] To begin 
with, we examine these images from a qualitative point 
of view. Starting from the case in which the Mn is in- 
serted in the second sublayer or deeper (Fig. [TT] (b)-(e)), 
the LDOS on the (110) surface displays a characteris- 
tic triangular shape with one vertex pointing down in 
the [001] direction. As the Mn is inserted into yet deeper 
layers, the triangular shape evolves into shape resembling 
a butterfly or bow-tie, with stronger upper wing. These 
features have been observed experimentally 7 8 and found 
in agreement with TB calculations 8 similar to ours. It 
is interesting to notice that for a Mn in bulk, calculated 
cross-section LDOS on (110) plane located n atomic lay- 
ers from the Mn shows a similar butterfly shape, but 
with a stronger lower wing. 8 Our calculations for bulk 
Mn shown in Figj4] displays a similar tendency - see 
in particular Fig. pT(d)-(e). This symmetry reversal of 
the bow-tie shape with respect to the [001] direction for 
a Mn close to the surface has been ascribed to the in- 
trinsic strain associated with the buckling relaxation. 8 
The change of the Mn-induced LDOS from a triangu- 
lar to a bow-tie shape as the impurity is inserted into 
deeper layers below the (110) surface occurs after sub- 
layer 5, where the magnetic anisotropy landscape starts 
to develop a quasi-easy plane. The Mn spin is now free 
to explore the easy-plane (110), with in more spectral 
weight in the [110] direction in the Mn layer below the 




FIG. 11: (Color online) The (110) surface LDOS as a function 
of the depth of Mn sublayer, (a)-(i) correspond to sublayers 1- 
9. Maximal intensity projections along [110] in the left column 
show how the LDOS in the ground state solution eventually 
detaches from the surface, suppressing the extension along 
the [110] direction. The maximal intensity projections show 
the relative distribution of the spectral weight in the cluster 
(thermometers apply to center and right column). The two 
right columns show the surface LDOS in the easy and hard 
direction (see Fig. |9j). The hard direction generally has a 
higher LDOS maximum and its pattern shows more spectral 
weight on the [001] side with increasing depth. The inset in 
(i) shows the surface LDOS when the Mn spin is pointing in 
the [001] direction, which is very close in energy to the easy 
direction. 
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surface. We now proceed with a more quantitative evalu- 
ation of Fig. Ill} The LDOS of the first sublayer (Fig. [IT] 



(a) ) is highly localized with a maximum spectral weight 
on the surface As of 50%. The Mn, one layer below, has 
only 8% spectral weight. The change in the LDOS with 
the Mn spin direction is minute and the resulting very 
low variations in orbital and spin character yield a low 
anisotropy energy. It is clear that the rising anisotropy in 
sublayer 2 and below is associated with a generally more 
spread out acceptor wavefunction. In sublayer 2 (Fig. [IT] 

(b) ), the hole is now much more spread out than at the 
surface and subsurface, and it shows that a large part of 
the extended spectral weight lodges at the surface. In 
sublayer 1 and 2 there is little change between the hard 
and easy direction, and in sublayer 3 and 4 the maximum 
spectral weight going from the hard to the easy direction, 
decreases by 14% and 33%. In sublayer 5 (Fig. [TT| (e)), 
where the maximum anisotropy energy is reached, the 
easy direction surface LDOS instead becomes more pro- 
nounced and the maximum LDOS decreases by 60% from 
the easy to the hard direction. There is clearly some- 
thing special about sublayer 5, where the surface has a 
high impact on the anisotropic extensions of the acceptor 
wavefunction. In addition to the large change in magni- 
tude between the hard and easy directions when the Mn 
is in this particular layer, the surface LDOS begins to 
show a qualitative change. In sublayer 6 (Fig. [IT] (f)), 
at which the high and low barriers are first interchanged, 
the wavefunction again shows a stronger surface LDOS 
in the hard direction. As we move further below the 
surface, the LDOS maximum decreases by around 40% 
as the Mn spin reorients from the hard to the easy di- 
rection. At the deepest level in sublayer 9 (Fig. 11 (i)) 



the quasi-easy (110) plane has formed, and the surface 
LDOS maximum abruptly drops 82% between the hard 
and easy direction. The general trend in acceptor wave- 
function character as layer depth is increased is that the 
surface LDOS is decreases and extends more along the 
[110] direction, running below the surface. 

Although the surface LDOS is biased on the [001] side 
of the Mn, the patterns are consistent with the fact that 
the wavefunction tends to extend along directions per- 
pendicular to the Mn spin. In the hard directions (where 
the Mn spin is in the [111] or the [110]), we therefore 
consistently see more spectral weight on the [001] side 
of the Mn relative the easy direction LDOS, and a more 
pronounced bow-tie like pattern. For the deeper levels 
where the low barrier has dropped significantly, thermal 
and quantum fluctuations, as well as Mn-Mn interactions 
in the sample can cause the Mn spin to fluctuate in the 
(110) plane, such that the bow-tie shape becomes more 
pronounced. The small inset in Fig. 11 (i) shows the 
LDOS at the surface when the Mn spin is in the [001] di- 
rection, which has the familiar slightly asymmetric bow- 
tie shape. This solution is very close in energy to the 
easy direction, only 1 meV higher. For deep impurities, 
the wavefunction is essentially dominated by the subsur- 
face extension along the [110], such that the Mn spin can 
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FIG. 12: (Color online) The LDOS of every atom in the clus- 
ter for the acceptor level for different Mn depths (marked 
0-4,9). The graph shows the 3200 values of LDOS (one for 
each atom) sorted in size and plotted on a logarithmic scale. 
The surface (marked 0) and the subsurface (1) both show a 
similar, highly localized distribution with a large population 
that has very low spectral weight in the range 10" 8 - 10" 7 . 
As the Mn is placed in deeper layers, the acceptor LDOS be- 
comes more extended with a majority population of atoms in 
the range 10" 4 - 10" 6 . 



move across the low barrier in the (110) plane. Com- 
paring with the fully periodic system, where the wave- 
function extends equally along the [110] and the [110] in 
the ground state, we see that the effect of the surface 
on the deep impurities is to reduce the extension along 
the surface normal, such that reversal paths open up in 
the single barrier for bulk (see Fig. [3|. It is also note- 
worthy that in sublayer 5, the low and the high barrier 
has not yet interchanged, but are comparable in energy. 
The LDOS at the surface when the Mn spin is point- 
ing in the direction of the lower barrier is larger than 
the LDOS in the easy direction, then following the same 
qualitative pattern as all the other depths. This indi- 
cates that the wavefunction is very sensitive to the level 
dynamics as the gap between the highest occupied and 
the acceptor level closes in. The odd behavior of sublayer 
5 is associated with a quasi-degeneracy between the high- 
est occupied and the acceptor level, as indicated by the 
large variations 5-50 meV of the gap between them (see 
Fig. 10). This sensitive situation, causes a large total 
amount of spectral weight to be shifted on and off the 
surface between the hard and easy direction, such that 
large variations in orbital and spin character occur, yield- 
ing the high anisotropy. 



In Fig. 12 the acceptor level LDOS values on all 3200 
atoms are sorted in size and then plotted on a logarith- 
mic scale. Note that what is shown is not a continu- 
ous curve illustrating a parameter dependence but 3200 
closely spaced points, one for each atom. The reason 
that the values have been sorted in size is that we are 
only interested in the possible magnitudes of the atomic 
LDOS here. There is no spatial information in this plot, 
but it shows how the LDOS magnitides are spread over 
the population of atoms. The sum of the LDOS values 
is normalized to unity. This plot demonstrates that an 
increased depth of the Mn is associated with a delocaliza- 
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nected with the variations in orbital and spin character. 
Fig. 13 shows the maximum and minimum orbital and 



1 0.020 
c 0.015 
o 0.010 
c 0.005 
Q-0.000 



FIG. 13: The evolution of the maximum and minimum orbital 
and spin character of the acceptor level on the unit sphere as 
a function of Mn depth. Orbital and spin character fluctua- 
tions are small in the highly localized surface and subsurface 
wavefunctions. They then increase to very large variations in 
sublayer 5, after which they level out. The spin down charac- 
ter follows a similar pattern, but drops sharply after sublayer 
4. 



tion of the acceptor wavefunction, which manifests itself 
as a large population of atoms with increased spectral 
weight. As the Mn depth increases a large population of 
atoms emerges with approximately 2-3 orders of magni- 
tude larger LDOS relative the highly localized surface set. 
The surface (marked 0) and the first sublayer (marked 1) 
exhibit similar distributions, with a localized signature 
and a larger spectral weight in the very narrow high end 
limit on the far right. For the surface and subsurface, the 
bulk of the spectral weight is located in close vicinity to 
the Mn, which results in a large population of atoms with 
a much lower LDOS in the range 1(T 8 - 1(T 7 . For sub- 
layer 3 and deeper, the majority of atoms show increased 
weight in the approximate range 10 -6 — 10 -4 . This large 
population represents atoms farther away from the Mn 
core region, which means that the acceptor wavefunction 
is becoming much more extended away from the surface. 

The spatial spread of the acceptor wavefunction is con- 



spin down character on the unit sphere of the Mn mag- 
netic moment directions as a function of Mn depth for 
the acceptor level. Layer (the surface) and layer 1 
(the subsurface) show little variation. The difference be- 
tween maximum and minimum orbital character rapidly 
increases from sublayer 3 and reach high values in sub- 
layer 5 and 6, where the anisotropy energy is also large. 
As we move further down, the orbital p-characters level 
out towards a 25% maximum and minimum difference. 



As can be seen in Fig. 10 the acceptor level comes very 
close to the highest occupied level in sublayer 5, lead- 
ing to a quasi-degeneracy between these two levels in the 
hard direction. This leads to a large anisotropy energy, as 
the quasi-degeneracy is lifted and the total energy low- 
ered, when the Mn spin is pointing in the easy direc- 
tion. Associated with this quasi-degeneracy, is a mixing 
of the two levels and a change in the nature of the accep- 
tor wavefunction. This change manifests itself as large 
fluctuations in orbital and spin character of the wave- 
function in sublayer 5, as can be seen in Fig. 13 The 



quasi-degeneracy between the acceptor and highest oc- 
cupied level in sublayer 6, is also related to the change in 
orbital and spin character responsible for the interchange 
of the high and low barriers as we go from sublayer 5 to 
6 (see Fig.P(e-f)). The spin-down character variation 
(Fig. [l3| (c)jincreases steeply in sublayer 2 and reaches a 
maximum in sublayer 4, after which the difference drops 
quickly towards 5%. The anisotropy is also quite high in 
sublayer 4, but the orbital character variations are not 
as large as in sublayers 5 and 6. However, looking at 
the variations in spin-down character, we see that there 
is a sharp maximum on sublayer 4. This points to the 
crucial impact of small spin character fluctuations of the 
acceptor level on the magnetic anisotropy energy. 

Fig. [12] and [13] show that the acceptor wavefunction 
delocalizes with depth, leading to larger variations in or- 
bital and spin character as the spread increases. Because 
of the surface, the variations in spin and orbital character 
are related to the geometry of the system that is responsi- 
ble for redistribution of the LDOS between the hard and 
easy directions. The total redistribution is very large for 
the maximal anisotropy layer, essentially depleting the 
surface LDOS in the hard direction to the advantage of 
the LDOS in the easy direction. 



IV. SUMMARY AND CONCLUSIONS 

We have undertaken a study of the properties of a sin- 
gle Mn embedded in a large 3200 atom GaAs matrix. 
Our model is based on a kinetic-exchange tight-binding 
Hamiltonian that accounts for the polarization of the As 
p-electrons nearest neighbors of a Mn spin via an off-site 
exchange term, and includes the local atomic spin-orbit 
interaction as well as the Coulomb field from the Mn 
ion. The relaxation of the (110) surface was taken into 
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account since it plays an essential role^El i n the (110) 
surface electronic structure. 

In agreement with previous calculations 6 we find that 
the acceptor wavefunction is generally highly dependent 
on the Mn spin direction. An exception occurs for the 
cases of surface and subsurface layer Mn spin locations, 
for which we find a highly localized acceptor wavefunc- 
tion that does not spatially redistribute as the Mn spin 
changes direction. We find, in particular, that a sin- 
gle Mn in the (110) surface produces a highly local- 
ized acceptor-level wavefunction and that the level lodges 
deep in the gap due to the loss of coordination at the sur- 
face. A deep acceptor level at the surface is also found in 
experiment!^, although uncertainly about band-bending 
effects has been responsible for some confusion as to its 
energetic position. The anisotropy of the surface state is 
relatively low, as expected for a deep-gap state. The ac- 
tual depth of this state in our model depends sensitively 
on a purely phenomenological parameter in our calcula- 
tion, the off-site Coulomb correction, which needs to be 
reduced relative to its bulk value in order to reproduce 
the experimentally observed position in the gap. 

The acceptor level spatial structure is probed exper- 
imentally via the LDOS contributions to the surface 
layer. As the Mn is placed in deeper subsurface lay- 
ers, the surface LDOS first displays a triangular pattern, 
which then evolves into an asymmetric butterfly or bow- 
tie shape when the Mn is located still further from the 
surface. This finding is in qualitative agreement with 
experiment!^ We find in addition that the anisotropy 
energy grows with depth, as the acceptor wavefunction 
becomes more extended. The magnetic anisotropy easy 
direction is strongly affected by the surface, and we find 
that its presence tends to favor [111] as an approximate 
easy direction. The observed easy axis is associated with 
a tilt of the wavefunction due to the up and down shift 
of the surface As and Ga atoms. This symmetry is re- 
sponsible for the triangular pattern on one side of the 
Mn and persists to deep layers, after which the Mn spin 
becomes more free to move across the surface normal. 
By comparing a three-dimensionally periodic calculation 
with the corresponding wide surface slab calculations, we 
see that for the deeper Mn, the effect of the surface is to 
open up additional reversal paths in the bulk Mn impu- 
rity's single-barrier magnetic anisotropy landscape. This 
is because at the deepest layers, the ground state solu- 
tion yields an acceptor wavefunction that is suppressed 
along the surface normal [110], relative to the fully pe- 
riodic system. The deep layer wavefunction extends in- 
stead most strongly along the [110] direction, such that 
a quasi-easy [110] plane is formed. A strong magneti- 
cally and spatially anisotropic state is found in sublayer 
5, where a quasi-degeneracy between the acceptor and 
highest occupied level leads to a large anisotropy energy 
and a large shift of the LDOS at the surface between hard 
and easy directions. 

Our predictions can in principle be checked experimen- 
tally by manipulating the Mn magnetic moment direction 



with an external magnetic field. As a function of the Mn 
depth, effects could be visible as early as the third sub- 
layer, where our calculations indicate that the surface- 
layer LDOS for magnetization along the [111] hard di- 
rection should be larger than the surface layer LDOS 
when the magnetization is along the easy direction [111]. 
In sublayer 3 and 4 the maximum LDOS decreases by 
14% and 33% from the hard to the easy direction. The 
acceptor level becomes increasingly shallow with increas- 
ing Mn layer depth and approaches the valence band in 
sublayer 5. We observe a strong signature in sublayer 
5. At this point the surface-layer LDOS is instead larger 
for magnetization along the easy direction, with a 60% 
decrease in maximum LDOS from the easy to the hard 
direction. It should be noted however, that this situation 
depends very sensitively on the local environment since 
this signature is associated with a quasi-degeneracy be- 
tween the highest occupied level and the acceptor level. 
For instance, increased Mn doping can affect the layer 
index at which this change occurs. In sublayer 6, the 
high and the low barriers have interchanged, with the 
high barrier now in the [110] direction parallel to the 
surface. At this point the easy direction is still approx- 
imately the [111], but it then moves towards [110] with 
further increased depth. In addition to a large 40% de- 
crease in the maximum LDOS from the hard to the easy 
direction, our calculations indicate that a more symmet- 
ric bow-tie shape appears in the hard direction for layer 
6 and deeper. 



The Mn spin-orientation can be influenced by thermal 
and quantum fluctuations as well as by external mag- 
netic fields. We will address the influence of quantum 
fluctuations in a subsequent publication, and comment 
here only on the interplay between external fields and 
thermal fluctuations. At temperatures higher than the 
anisotropy energy the Mn spin orientation will be ran- 
domized by interactions with its thermal bath. The 
measured surface-layer LDOS should then correspond 
to an average of the results obtained here for partic- 
ular orientations. At low-temperatures thermal fluctu- 
ations become unimportant and the magnetization ori- 
entation will depend on a competition between Zeeman 
coupling and magnetic anisotropy. As we have explained, 
the surface-layer LDOS pattern provides an indirect fin- 
gerprint of the magnetization orientation. Our exper- 
imental predictions for the depth dependence of mag- 
netic anisotropy can be tested by identifying the Zeeman- 
coupling strength required to change the surface-layer 
LDOS. For example we predict that for Mn in deeper lay- 
ers there are high magnetic barriers of the order 10 meV, 
which implies that large magnetic fields are required to 
reorient the Mn spin to the hard directions. These high 
barriers for a single Mn suggests that it might be possible 
to engineer few-atom substitutional Mn impurity clusters 
near the GaAs surface which act like nanomagnets with 
attractively large magnetic blocking temperatures. 
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